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THE NON-POLYNOMIAL CONSERVATION LAWS AND INTEGRABILITY 
ANALYSIS OF GENERALIZED RIEMANN TYPE HYDRODYNAMICAL 

EQUATIONS 

ZIEMOWIT POPOWICZ AND ANATOLIY K. PRYKARPATSKY 

Abstract. Based on the gradient-holonomic algorithm we analyze the integrability property of 
the generalized hydrodynamical Riemann type equation D^u = for arbitrary A'^ S Z+. The 
infinite hierarchies of polynomial and non-polynomial conservation laws, both dispersive and 
dispersionless are constructed. Special attention is paid to the cases iV = 2, 3 and N = 4, 
for which the conservation laws, Lax type representations and bi-Hamiltonian structures are 
analyzed in detail. We also show that the case A'^ = 2 is equivalent to a generalized Hunter- 
Saxton dynamical system, whose integrability follows from the results obtained. As a byproduct 
of our analysis we demonstrate a new set of non-polynomial conservation laws for the related 
Hunter-Saxton equation. 



1. Introduction 



CN , Nonlinear hydrodynamic equations are of constant interest still from classical works by B. Rie- 

^T ' mann, who had extensively studied them in general three-dimensional case, having paid special 

attention to their one-dimensional spatial reduction, for which he devised the generalized method 
of characteristics and Riemann invariants. These methods appeared to be very effective [U H] [18] 
in investigating many types of nonlinear spatially one-dimensional systems of hydrodynamical type 
and, in particular, the characteristics method in the form of a " reciprocal" transformation of vari- 
ables has been used recently in studying a so called Gurevich-Zybin system [2 [3] in |8] and a 
Whitham type system in [18j|9]. Moreover, this method was further effectively applied to studying 
solutions to a generalized |10j (owing to D. Holm and M. Pavlov) Riemann type hydrodynamical 



system 



H 

P. ; (1.1) D^u = 0, Dt := d/dt + ud/dx, N e Z+, 

where dx/dt ^ u E C°° (M; M) is the corresponding characteristic flow velocity along the real axis 
M. 

We will consider, for convenience, the hydrodynamical equation (jl.ll) on the 27r-periodic space of 
functions A^o '■= C°°(M/27rZ;R), which can be, obviously, equivalently rewritten as the following 
nonlinear dynamical system in the augmented functional manifold M := C°°(R/27rZ;M^) on the 
vector u:= (u^ := u, u^^^ := DM°\ m(2) := dM^\ ..., u(^-i) := Dtu^^-^^y G M : 



4^) = z.(2)-J°)42)^ 



(1.2) 






The dynamical system (jl.2p possesses a very interesting and important property: the partial 
flows of the velocity components u'^-^\j ~ 0,N —1, are realized along the axis M with the same 
characteristic velocity dx/dt = m^°). This is exactly the case, deeply studied by Riemann (see, for 
example, [T]) when one can introduce the so-called "Riemann invariants" making it possible to 
obtain a suitable separation of dependent variables important for the integration. Really, we can 



(n\ 
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recurrent set of differential equations in full differentials; 

(1-3) 

du(^-i) = 0, 
whose solution is easily found by means of simple integration in the parametric form as 



(N-l) 






2! 



3! 2! ^ ' 



(1.4) 



u(i) 



^tA^-2 ^ p 



u 



(0) _ 






Y^ -jl3N^2~jiz), 



Tv + E-y/^^-i-^W' 



(A^-i)! ^J- 



where l3j e C°°(M;R), j = 1, A^- 1, and I3n G C°° (M/27rZ; K) are arbitrary smooth functions, 
depending on a suitable first integral z S C°° (K^ ; R) of the system (|1.3I) . The presented above 
result coincides in some part with that obtained before by M. Pavlov in ll4j, and can be used for 
constructing special solutions in analytical form to the generalized Riemann type hydrodynamical 
equation (|l.ip . 

As it was stated before in [161 IS H HSl HH [10] the Riemann type hydrodynamical system (|1.2I) at 
N — 2 and A^ = 3 possesses additional very interesting properties, being a completely integrable bi- 
Hamiltonian system. In particular, it possesses infinite hierarchies of dispersionless and dispersive 
conservation laws, which can have an important hydrodynamical interpretation and may be used 
for constructing a wide class of other special quasi-periodic and solitonic solutions. 

In spite of the exact integrability of dynamical system (jl.ip by means of the classical charac- 
teristics method, the solutions obtained this way are, to the regret, of very vague usefulness, as 
they are given in the entangled and involved form not fitting for studying solutions belonging to 
some specially assigned classes of functions, for instance, fast-decreasing, quasi-periodic and etc. 
Thereby, additional studying of the mathematical structures associated with dynamical system 
p.ip by means of modern symplectic theory techniques is as much as could needed, and what is a 
topic of our present investigation. 

The next section below is devoted to the Hamiltonian analysis of the hydrodynamical system 
(|1.2p at A^ = 2, A^ = 3 and A^ = 4, as well as to the description of their new hierarchies of 
conservation laws, the related co-symplectic structures and Lax type representations. 



2. The generalized Riemann type hydrodynamical equation at N=2: conservation 

LAWS, BI-HAMILTONIAN STRUCTURE AND LAX TYPE REPRESENTATION 

Consider the generalized Riemann type hydrodynamical equation (jl.ip at A^ = 2 : 
(2.1) D^u = 0, 

where Dt = d/dt + ud / dx, which is equivalent to the following dynamical system: 

2.2 > := K\u,v\, 
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where K : M ^ T{M.) is a related vector field on the 27r-periodic smooth nonsingular functional 
phase space M := {{u,vY e C°°{R/2TrZ;R'^) : ul - 2vx ^ 0,x £ M}. As we are interested first in 
the conservation laws for the system (|2.2p . the following proposition holds. 

Proposition 2.1. Let H(X) := Jj^ h(x; X)dx G D{M) be an almost everywhere smooth functional 
on the manifold Ai, depending parametrically on A G C, and whose density satisfies the differential 
condition 

(2.3) ht = X{uh), 

for all t (zM. and A G C on the solution set of dynamical system 112. 2\) . Then the following iterative 
differential relationship 

(2.4) if/h)t - Xiuf/h), 

holds, if a smooth function f G C°°(]R;R) (parametrically depending on A G C) satisfies for all 
t G M the linear equation 

(2.5) ft^2Xu.J + Xuf,. 
Proof We have from ^^^-l^^b^ that 



(2.6) {f/h)t = ft/h- fht/h" ^ ftlh-Xfu^/h-Xfuh^/h^ = 

= ft/h + Xfu{l/h),^XuJ/h = 

= X{uf),/h + Xuf{l/h),^X{uf/h),, 

proving the proposition. D 

The obvious generalization of the previous proposition is read as follows. 
Proposition 2.2. // a smooth function h G C°°(K;]R) satisfies the relationship 

(2.7) ht — kuxh + uhx, 
where fc G M, then 

/.27r 

(2.8) H= h^/'^dx 

Jo 

is a conservation law for the Riemann type hydrodynamical system 112. 2\) . 

Remark 2.3. Let h G C°°(M;R) satisfy the differential relationship ht — {hu)^, then / = /i^ is a 
solution to equation p.4p . 



Remark 2.4. If functions hj G C°°(M;M),j G Z+, satisfy the relationships hj^t — ^{hju)xTJ G ^+j 
A G C, then the functionals 

(2.9) F(,,)= Y. ki:-^' Th^^h!^-'^^ 

with kn G M, n G Z+, z,j G Z_(_, being arbitrary constants, are conserved quantities to equation 
()2.2p . This formula, in particular, makes it possible to construct an infinite hierarchy of non- 

rirtKrnnmicil rrwiaorTrtirl mianfi-l-ifiG fnr -j-lnti R i i=>m a n n -j-^r-nti VnT-rlmrl^rnQmir'al G-\7-G-ftim (\0 9I\ 
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Example 2.5. The following non-polynomial functional 



Hi^^ 



Jo Jo 

h\^^ = f ^ villi - 2vx)dx, 
Jo 

i?g^ = / {kiu{uxxVx - UxVxx) + kiVxxV + k2{ul.v - 2vl.)) dx, 



(2.10) 



(1) _ z-^- 1 

^Q / y^xx'^xxx ^xxx'^xx ) (xX, 

Jo 

-ffg" = / {Ux{UxxVx - UxVxx) + VxxVx)Y , 

Jo 

2tt 



J^IO — / \'^'^xx[^x'^xx '^xx'^xj ^xx) 



_ _ 2 ^(6 

y^^xxy^x^xx ^xx^xj ^ 


are conservation laws for the Riemann type dynamical system p.2p . 

Quite different conservation laws have been obtained in [71 [8] using the recursion operator 
technique. The corresponding recursion operator proves to generate no new conservation law, if 
one applies it to the non-polynomial conservations laws ()2.10p . 

We also notice that dynamical system (|2.2p . as it was shown before in [5J[T3], can be transformed 
via the substitution 

(2.11) y^^d'\ul + rj^) 

into the generalized two-component Hunter - Saxton equation: 

(2-12) Ux,t = --ul-uuxx + -^r]'^, 

Vt = -iuv)x- 

This equation allows the simple reduction to the Hunter - Saxton dynamical system [5l [151 [13] at 
77 = 0: 

(2.13) Uxt^--;:ul-uuxx- 

The non-polynomial conservation laws (j2.10p . upon rewriting with respect to the substitution 
(12.111) , give rise to the related non-polynomial conservations laws for the generalized two-component 
Hunter - Saxton dynamical system ()2.12[) . Moreover, if we further apply the reduction rj = 0, we 
obtain, respectively, new non-polynomial conservation laws for the Hunter - Saxton dynamical 
system (|2.13p , supplementing those found before in [TSl US] ■ 

Example 2.6. The following functionals 

^0 1/1^ u(^) r^" f 2\^J ud) f^"" UxxxU + 2UxxUx , 

(2.14) H^-" = / [uxxK) dx, i?9 = / -^= dx, 

Jo Jo \/Uxx 

1 r^'" 1 

Hs^ = / [u^xUxid^^ul) -Uxxulul^dx 
Jo 

are the conservation laws for the Hunter-Saxton dynamical system (|2.13p . 

All of these and many others non-polynomial conservation laws can be easily obtained using 

Ti :j.: /TTTTli t:i i_ ii i r_ j.: l^ 
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Jo Jo 

(2.15) iJ(2) == / y/u^dx, iJ(3) = / y/u^^{d-^ul)dx, 

Jo Jo 

F(4) = f \id-'ul){uU.,ul,-uUd-'ul))]hx 

Jo 
are also conservation laws for the Hunter-Saxton dynamical system (j2.13p . where m ^ —An and 
n, m € Z. 

Now we proceed to analyzing the Hamiltonian properties of the dynamical system (j2.2p . for 
which we will search for solutions to the determining [THl HI] Nother equation 

(2.16) Lk^ = ^t- ^K''* - K'd = 0. 

where Lk denotes the corresponding Lie derivative on M. subject to the vector field K : M. ^ 
T{M), K' : T{M) -^T{M) is its Frechet derivative, K'* : T*{M) ->T*{M) is its conjugation 
with respect to the standard bilinear form (•, •) on T*(A^)xr(A^), and d : T*{M) -^ T{M) is a 
suitable implectic operator on M , with respect to which the following Hamiltonian representation 

(2.17) K^-d grad Ha 

for some smooth functional H^ G D{Ai) holds. To show this, it is enough to find, for instance by 
means of the small parameter method Ll8j.20^, a non-symmetric {ip' ^ ip''*) solution ip g T*{Ai) 
to the following Lie-Lax equation: 

(2.18) i^t+K''*^^ grad C 

for some suitably chosen smooth functional C £ 'D{M). As a result of easy calculations one obtains 
that 

1 r^"" 

(2.19) ip = (v,0)T, L = - v^dx. 

2 Jo 

Making use of p.lSp jointly with the classical Legendrian relationship 

(2.20) H^ := (V^, K) - C 

for the suitable Hamiltonian function, one easily obtains the corresponding symplectic structure 

(2.21) ,?-!:= V'-V'''*=(_°1 J) 
and the non-singular Hamilton function 

1 r^'^ 

(2.22) Ha- - [v"^ +v^u^)dx. 

2 Jo 

Since the operator (j2.2ip is nonsingular, we obtain the corresponding implectic operator 

(2.23) ..(; -1 

necessarily satisfying the Nother equation (|2.16p . 

It is worth here to observe that the determining Lie-Lax equation p.lSp possesses still other 
solution 

(2.24) ^^(_iL,„^), £= / uv^dx, 

2 2vx 4 Jo 

giving rise, owing to expressions (I2.2ip and (12.201) . to the new co-implectic (singular "symplectic") 
structure 



(2.25) 77-1 := 7/.' - V''* - ' "" '^^xv.. 



on the manifold A4, subject to which the Hamiltonian functional equals 

1 f^- 
(2.26) H„ := - / (ut.v — Vxu)dx, 
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supplying the second Hamiltonian representation 

(2.27) K ^ -T] grad H^ 

of the Riemann type hydrodynamical system p.2|) . The co-implectic structure (|2.25p is singular, 
since ?)~^(wa;, ^a;)^ = 0, nonetheless one can calculate its inverse expression 

Moreover, the corresponding implectic structure 77 : T*{Ai) -^ T*{Ai) satisfies the determining 
Nother equation 

(2.29) LkV = vt- VK''* - K'rj = 0, 

whose solutions can also be obtained by means of the small parameter method [T51 [TH]- We 
remark also that, owing to the general symplectic theory results ]TE[ UM [201 HI] for nonlinear 
dynamical systems on smooth functional manifolds, operator (|2.25p defines on the manifold A4 
a closed functional differential two-form. Thereby it is a priori co-implectic (in general, singular 
symplectic), satisfying on M the standard Jacobi brackets condition. 

As a result, the second implectic operator (|2.28|) . being compatible [T8j [21] with the implectic 
operator ()2.23p , gives rise to a new infinite hierarchy of polynomial conservation laws 



(2.30) 



I dX< {d^'^rjygrad H^[uX},u> 
Jo 



for all n G Z+. Having defined the recursion operator A :— ■& ^77, one also finds from p.30|) . (|2.16|) 
and (|2.29p that the following Lax type relationship 

(2.31) LKA^At-[A,K''*]^0 

holds. If to construct now the asymptotical expansion 1^9(2;; A) ~ ^ X^^^^ grad 7j_i[u,w] as 
A -^ 00, it is easy to obtain from (j2.30p that the gradient like relationship 

(2.32) X'^^^{x;X^'nip{x;X) 

holds. The latter relationship, making use of the implectic operators (|2.23l) and (I2.28p . can be 
represented in the following two factorized forms: 

M^-A) \^{ -4A3/2 + 2A«,/| \ / -2A(/i/2). 
<^2(a;;A) i [ -4Xy,h~2Xu.,fi J I -(/|), 



(2.33) ^(.;A):= ^^.^ = -^ , ,_ ,,^, ,, , - , ,,,, 



where a vector / S C°°(R^; C^) lies in an associated to manifold M vector bundle £(A^; E^), whose 
fibers are isomorphic to the complex Euclidean vector space E^ . Take now into account [18l [20] 
that the Lie-Lax equation 

(2.34) LKifix; A) = dip{x; X)/dt + K'^*ip{x; A) = 

can be transformed equivalently for all a;, i G M and A G C into the following evolution system: 

(2.35) Dtf =( ^. ^"^ ^V, Dt= d/dt + ud/dx. 

\ -L 7/2; J 

The equation (|2.35p . owing to the relationship (|2.32p and the obvious identity 

(2.36) Dtf^ + u,f,^{Dtfy„ 
can be further split into the adjoint to (|2.35p system 





(2.37) Dtf - <z(A)/, qiX) ■■= ^ _^ Q 
where a vector / G C°°(M^; C'^) satisfies the following linear equation 

(2.38) L^£\u^v:X]f, i\u,v;X] -.^ ( '}}t'' 7''"^ 
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compatible with ()2.37|) . Moreover, as a result of l\2.37\i and (|1.4p . the general solution to (|2.38|) 
allows the following functional representation: 

(2.39) fiix,t) ^ gi{u~tv,x-tu + vt^/2), 

f2{x,t) = -t\gi{u-tv,x-tu + vt'^/2) + 
+92{u — tv,x ~ tu + wt^/2), 

where gj G C°° {M.'^ ; C) , j = 1,2, are arbitrary smooth complex valued functions. Now combining 
together the obtained relationships ()2.37|) and (|2.38p . we can formulate the following proposition. 

Proposition 2.7. The Riemann type hydrodynamical system il.l\) is equivalent to a completely 
integrable hi-Ham.iltonian flow on the functional manifold A4 , allowing the Lax type representation 



f, = l[u,v-,\\f, !,=p(t)f, f{l):=-ul[u.v-\\+q(\), 

Z ) ■ '<*) - 

XUxU V,j;U 



(2.40) e[u,v;X]:=( ^^,^ ^^A , q{\) := ( \ [] 



^'^^^ \ -A-2A^w -Xu^u 

where / £ C°°(M^; C^) and A e C is an arbitrary spectral parameter. 

Remark 2.8. It is worth to mention here that equation (12.37^ is equivalent on the solution set of 
the Riemann type hydrodynamical system (12.21) to the alone equation 

(2.41) Dlf2 = ^^ A/i = 0, A/2 = -A/i, 

where vector / G C°°(M^; C^) satisfies for all A G C the compatibility condition (12.381) and whose 
general solution is represented in the functional form (|2.39p . 

Concerning the set of conservation laws {TJq ' \HI ' '}, constructed above, they can be ex- 
tended to an infinite hierarchy {iJ- € D{M) : j e Z+ U { — 1}}, where 



(1/2) ._ 



(2.42) H^'''> :^ / a2j+i[u,v]dx 



2tt 







and the affine generating function <T(a;;A) :— ^ln/2(x;A) ~ X]i=z u{-i> '-"^it"'^]'^ ^ as A — J- oo 
satisfies the following functional equation: 

(2.43) {a- Xu^)^+a^ + X^{2v^-ul) = 0. 



271- 



In addition, the gradient functional (p{x; A) :— grad ^{x; A) G T*{M), where 7(A) := Jg a{x; X)dx 
satisfies for all A G C the gradient relationship (|2.32p . 



3. The generalized Riemann type hydrodynamical equation at N=3: conservation 

LAWS, BI-HAMILTONIAN STRUCTURE AND LAX TYPE REPRESENTATION 

Here we proceed to analyzing conservation laws and bi-Hamiltonian structure of the generalized 
Riemann type equation (jl.ip at A^ = 3 : 

Ut — V — UUx 

(3.1) vt^z~uvx ):^K[u,v,z\ 

Zt = -UZx 

where K : M. ^^ T{M.) is a suitable vector field on the periodic functional manifold M := 
C°° (R/2ttZ;M.^) and t G M is an evolution parameter. The system p.ip proves also to possess 
infinite hierarchies of polynomial conservation laws, being suspicious for complete and Lax type 
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Namely, the following polynomial functionals are conserved with respect to the flow (|3.1 

(3.2) ijW : = f\xz"{vu,-v,u^^z), 

Jo " + 1 

p2TT 

= / dx[—7vxV u + z{6zu + 2vxU — 3v — 4vuux)], 
Jo 



ij(4) 



p27T plTT 

ij(5) : =/ dx{z^ux~2zvvx), H^^'> -.^ dx{z^v^ + Sz^v^u + zJ), 

Jo Jo 

ij(^) : =/ dx{zxv'^ + iz^vux-iz'^), 
Jo 

n2-K 

H^ ' : = / dxz{6z u + 3zVxU — 3zv — AzvUx — 2vxV u + 2v Ux), 
Jo 

/■27r 

ij(^) : =/ dx[1001vxv\+{1092z'^u^ + 36AzVxU^ - 
Jo 

-1092zv'^u - 728zvuxu'^ - MAyxV^u"^ + 273w'* + 728v^Uxu]), 

p2ir p27T 

Hi^'> : = dxzxvz'\ H^^') := dxzx{v^ - 2zuY , 

Jo Jo 

where n G Z+. In particular, as n = 1, 2, ..., from (|3.2p one obtains that 

/•2lT r2iT 

(3.3) 



r(2) . 




= / dxZxV, H\ :— 


/ dxzxzv,..., 




Jo 


Jo 




.27r 




r(3) . 
1 


= / dxZx(v — 2uz), 
Jo 

<.2ir 








r(3) . 

2 


= / dxzxiv"^ + 4z^u^ 
Jo 


-4zu^u),..., 



and so on. 

Making use of the iterative property, similar to that, formulated above in Proposition 12. H one 
can construct the following hierarchy of non-polynomial dispersive and dispersionless conservation 

laws: 

<.27r 
11^ =^ I uXy ZUxx^x-^x ~r ^xx'^x ' ^^x^xx 

Jo 



^x^xx'^x ~r OVxx^x "^^x^xx ) 5 
(.27r 
1I2 / CLXy Vxx'^x I ^x^xx) 1 

p2-K 
-^^3 / ^XyVxx'^x ^x'^xx ^xx) 7 

Jo 

/'27r 

iyf''^^ - / dx[~2vuxzx+vl + z{-uxvx + 3zx)]^/^, 
Jo 

<.27r 

(3.4) H^^^^'' = / dx{8ulzx~3ulvl-l8uxVxZx + Qvl + 9zx)^^'^, 



Jo 

p2-K 
rl-^ — / dXy Zllxxx^x^x ~r ^xxx^x ' ^xx^^ ^^xx^x^xx ~ 

Jo 

'jUxx^xx^x 'y ZUxZxxx '^x'^xxx'^x I 

-\-oV,xVxx ' ^^XXX^X '^'^X^XXX ) J 

p2ir 
Hi = / dx[kiu{~VxxZx + VxZxx) + kiv{UxxZx - UxZxx) + 



'3 

/O 



^z{k2UxxVx - k2UxVxx + kiZxx + fe^aja;) + 



THE NON-POLYNOMIAL CONSERVATION LAWS AND INTEGRABILITY ANALYSIS OF GENERALIZED RIEMANN TYPE HYDRODYNAMICA 

where kj G M, j' = 1, 3, are arbitrary real numbers. Below we will attempt to generalize the crucial 
relationship (|2.37p from Section 2 on the case of the Riemann type hydrodynamical system (I3.ip . 
Namely, we will assume, based on the Remark (j2.3|) . that there exists its following linearization: 

(3.5) D^MX) = 0, 

modeling the starting generalized Riemann type hydrodynamical equation (jl.ip at TV = 3, and 
where /3(A) g C°°(M^;C) for all values of the parameter A G C. The scalar equation p.Sp can be 
easily rewritten as the system of three linear equations 

(3.6) A/i = 0, A/2 = Mi/i, A/3 = M2/2 

where we have defined a vector / := (/i, /2, /2)^ G C°°(M'^; C'^) and naturally introduced constant 
numbers /ij := /ij(A) G C, j ~ 1,2. It is easy to observe now that, owing to the former result 
()1.4p , the system of equations p.6p allows the following solution representation: 



fi{x,t) = giiu-tv + zt^/2,v- zt,x-tu + vt^/2- zt^/6), 
f2{x,t) = tfiigi{u-tv + zt'^/2,v- zt,x~tu + vt'^/2- zt^/6) + 

(3.7) +52(w -tv + zf/2, V- zt,x-tu + vt^/2 - zt^/6), 

f^{x,t) = ^i^2—gi{u~'tv + zt^/2,v-zt,x-tu + vt^/2-zt^/Q) + 

+tii292{u -tv + zt'^/2, V- zt,x-tu + vt^/2 - zt^/6) + 
+gz{u -tv + zt'^/2, V- zt,x-tu + vt'^/2 - zt^/6), 

where gj G C°°(R'^; C), j = 1, 3, are arbitrary smooth complex valued functions. The system p. 6 
transforms into the equivalent vector equation 

/ 

(3.8) Dtf^qifi)/, g(A):= a*i(A) 

V /i2(A) 

which should be compatible both with a suitably chosen equation for derivative 

(3.9) U=e[u,v,z;X]f 

with some matrix £[u,v, z; X] G SL{3;C), defined on the functional manifold M, and with the 
Lie-Lax equation ()2.34p . rewritten as the following system of equations 

(3.10) Dt^ = 

where the vector ip := (p{x; A) G T*(A4) is considered as the one factorized by means of a solution 
/ G C°°(M2;C3) to ([3^ . satisfying the identity (f236l) . Namely, it is assumed that the following 
quadratic trace-relationship 

(3.11) ip^tri^X;u,v,z]f(g>P) 

holds for some vector valued matrix ^[X;u,v,z] G E'^ ig) End E^, defined on the manifold Ai, 
where "(g)" means the standard tensor product of vectors from the Euclidean space E''. Making 
now use of the determining expressions (j2.36p . (j3.1ip and p.Sp . one can find by means of some 
slightly cumbersome but tedious calculations that /ii(A) ~ A, fJ-2iX) = 1, A G C, and the matrix 
representation of the derivative p.9p 

(Xux 
3A2 
QX'^r[u, V, z] 

compatible with determining equation (|3.10p . where a smooth mapping r : A4 ^> M. satisfies the 
differential relationship 






Vx 


Zx 


\ 






-1 


— Ux 





u, 


Dt 


= d/dt + ud/dx 





-1 


-Ux 


/ 
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The latter possesses a wide set TZ of different solutions amongst which there are the following: 
(3.14) r € 7^ := {[{xv - uy2)/zU {v^ - ul/6)z-\ "'/\~ "-"-_+f ^/^ 

{VxV /G^UxV z/2 + UZx{uZ — V )/6 + vz )z^ }. 

Note here that only the third element from the set p.l4p allows the reduction z = to the case 
N — 2. Thus, the resulting Lax type representation for the Riemann type dynamical system p.ip 
ensues in the form: 

fx^£[u,v,z;X\f, ft^p{£)f, p{i) ■.^~ue[u,v,z;X\+q{X), 



XUx 

i[u,v,z;X]= \ 3A2 -2Xux Xvx I, ^(A) 





(3.15) 



— XuUx UVx —UZx 

p(£) — — 3mA^ + A 2XuUx —XuVx 

-6A^r[w, u, z]u 1 + 3uA —XuUx 

where / € C°°(R^; C^) and A G C is a spectral parameter. 

The next problem, which is of great interest, consists in proving that the generalized hydro- 
dynamical system (|3.ip is a completely integrable bi-Hamiltonian flow on the periodic functional 
manifold A4, as it was proved above for the system (j2.2p . 

That dynamical system p.ip is bi-Hamiltonian that follows easily as a simple corollary from 
the fact that it possesses the Lax type representation (|3.15l) and from the general Lie-algebraic 
integrability theory [171 IHl HI] ■ Taking into account that dynamical system (|3.ip possesses many 
(at least 4) Lax type representations, one derives that it possesses many (at least 4) different pairs 
of compatible co-symplectic structures, every of which generates its own infinite hierarchy of com- 
muting to each other conservation laws. Moreover, the involution of conservation laws belonging 
to different hierarchies fails owing to their non-compatibility. As finding of these structures is 
adjoint with cumbersome enough analytical calculations, we present below only a one pair of re- 
lated co-symplectic structures, making use of the standard properties of determining them Lie-Lax 
equation (|2.18|) . 

To tackle with the related task of retrieving the Hamiltonian structure of the dynamical system 
p.l[) . it is enough, as in Section 2, to construct [ISl [5D] exact non-symmetric solutions to the 
Lie-Lax equation 

(3.16) i;t + K''*il; = grade, V'VV'''*, 

for some functional C € D{Ai), where ip e T*{A4) is, in general, a quasi-local vector, such that 
the system p.ip allows the following Hamiltonian representation: 



K[u,v,z] = -r] grad H[u,v,z], 

As a test solution to (|3.16p one can take the one 

ip = {ux/2,0,-z^^ul/2 + z^^VxY , ^=0 {2z + vux)dx, 

^ Jo 



which gives rise to the following co-implectic operator: 
(3.18) 77-^ := V - V''* = 



f d -duxz-^ \ 

dzx 



-{ulzx'^d + dulz,^)- 

V ""^^ "" ^^'^ -{VxZ-^d + dVxZ-^) ) 



-UxZ^ ^d Zxd 2 



This expression is not strictly invertible, as its kernel possesses the translation vector field d/dx : 
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Nonetheless, upon formal inverting the operator expression p.l8|) . we obtain by means of simple 
enough, but slightly cumbersome, direct calculations, that the Hamiltonian function equals 



(3.19) 



Hr. 



dx(uj.v 



and the implectic ?7-operator looks as 



(3.20) 



a-i u^a-i 

Zrd-^ 



The same way, representing the Hamiltonian function (I3.19P in the scalar form 



(3.21) 



Hri ^ {■lp,{Uj:,Vx,Zxy), ■,p^-(-v,U + 



--d~'V~z)\ 



one can construct a second implectic (co-symplectic) operator d : T*{jv{) — > T(A^), looking up to 
0{fj.^) terms, as follows: 

(3.22) 



/ 



f((^9 + .W9) 



f(^^ + ^^)+ 

+f (wWa + au(i)) 



f(9(^ + 5.W) \ 

/i(az(i) + z(i)a) J 



where we put, by definition, "d^^ := {ip' — ^''*), u :— ijlu^^\v := iiv^^\z :— jiz^^' as /i — > 0, 
and whose exact form needs some additional simple but cumbersome calculations, which will be 
presented in a work under preparation. 

The operator p.22p satisfies the Hamiltonian vector field condition: 

(3.23) {ux.v^.z^y ^ -■& grad Hj^, 

following easily from p.2ip . 

The results obtained above can be formulated as the following proposition. 

Proposition 3.1. The Riemann type hydrodynamical system hl.l]) at N = 3 is equivalent to 
a completely integrahle hi- Hamiltonian flow on the functional manifold Ai, allowing the Lax type 
representation I13.15\) and the compatible pair of co-symplectic structures 113. 20\) and 113. 22\) . 



The infinite hierarchy of conservation laws like p.4p and related recurrent relationships can be 
regularly reconstructed, if to compute the asymptotical solutions to the following Lie-Lax equation: 

Lj^(f = 1^^. + K''*ip ^0, 

(3.24) ^ ~ a{x;X)exp{X^T + d-^a{x;X)}, 

where, by definition, a{x;X) ~ X^/ez o,j[u,v,z]X^^ , a{x;X) ~ J^jei, ui-n^ii'^^^^'^]^"^ as A — >■ 
oo, and 



(3.25) 



rfr 



{u, V, zy 



(1/3) 



-3ry grad TJg [u, v, z] = 

-{ulh-%+v-'{vlh-% 
-VxU^Hulh ^)^ + z^^{zlh 2)^ } ■.^K[u,v,z], 

— ZxU^ (Z^fl )x 

h{u, V, z\dx, h[u, V, z] = {vx.xUx - u^^v^ - z^^yl'^. 



H. 



is a Hamiltonian vector field on the functional manifold AA with respect to a suitable evolution 
parameter r e M. Since the vector fields p.25p and p.ip are commuting to each other on the 



whole manifold AA , the functionals H 



(1/3) 



1 J_l_ J_l_ 



Jq dj[u^v^z\dx^ J G Z+ U { — 1}, will be functionally 



-0(^2), 
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The Lax type integrability of the Riemann type hydrodynamical equation (jl.ip at A^ = 2 and 
A^ = 3, stated above, allows one to speculate that it is also integrable for arbitrary N E Z+. 

Concerning the evident difference between analytical properties of the cases N = 2 and A^ = 3, 
we can easily observe that it is related with structures of the corresponding Lax type operators 
(j2.38p and (|3.15p : in the first case the corresponding r-equation (|3.13p is trivial (that is empty), 
but in the second case it is already nontrivial, allowing many different solutions. This situation 
generalizes, as we will see below, to the case N > 4, thereby explaining the appearing diversity of 
related Lax type representations. 

To support this hypothesis we will prove below that also at TV = 4 it is equivalent to a Lax 
type integrable bi-Hamiltonian dynamical system on the suitable smooth 27r-periodic functional 
manifold M := C°°(M/27rZ;M'*), possesses infinite hierarchies of polynomial dispersionless and 
dispersive non-polynomial conservation laws. 



4. The generalized Riemann type hydrodynamical equation at N=4: conservation 
LAWS, bi-Hamiltonian structure and Lax type representation 

The Riemann type hydrodynamical equation (11.11) at A^ = 4 is equivalent to the nonlinear 
dynamical system 



(4.1) 



Ut = V - 


- UUr^ 


Vt = w - 


- UVa 


Wt=Z- 


- UW 



:= K[u,v,w,z], 



zt 



where K : M ^ T{Ai) is a suitable vector field on the smooth 27r-periodic functional manifold 
A4 := C°°(M/27rZ; M^). To state its Hamiltonian structure, we need to find an exact non-symmetric 
functional solution ip E T*{A4) to the Lie-Lax equation p.l6|) : 



(4.2) il^t+K' ip = grad C 
for some smooth functional £ E D{M), where 

/ -~du 1 0^ 

(4.3) k' = 



V 



-Z:r 



1 

-ud 







1 

-ud 





1 
-ud ] 



K = 



/ ud 
1 


V 



-Zx \ 



du J 

T{M) and its conjugate. The 



du 
1 






du 
1 



are, respectively, the Frechet derivative of the mapping K : Ai 

small parameter method [18^, applied to equation (|4.2p . gives rise to the following its exact solution 

„,2 ,, , 



(4.4) 



^ = 



ZZx 



{zux — vwxli)dx. 



-y 



As a result, we obtain right away from (|4.2|) that dynamical system (|4.ip is a Hamiltonian system 
on the functional manifold M., that is 



(4.5) K 

where the Hamiltonian functional equals 



-'d grad iJ, 



(4.6) 



H := (Va K)- £ 



2-K 



and the co-implectic operator equals 

/ 



(4.7) 



^- 



V-' - i)'' 



\ 









~d 





—ud 





-9 








w^Q 


_^a 


u^Q 



{uZx — VWx)dx 



-d^ 

l[z-^{vl'-2uxWx)d+ 



/ 



-{-divl ~ 2uxWx)z^^] 
The latter is degenerate: the relationship ■d~^{ux,Vx,Wx, z^y = holds exactly on the whole 



;r_i J A .rf T__ 



J ] 1 __i i^j ] i_ 
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To state the Lax type integrability of Hamiltonian system (|4.ip we will apply to it, as in Section 
3 above, the standard gradient- holonomic scheme of |18[[20] and find the following its linearization: 

(4.8) 



DthW - 0, 



where /4(A) £ C°°(M^; C) for all A £ C. Having rewritten (14.81) in the form of the linear system 



(4.9) 



A/ = g(A)/, q{X) 



^0000^ 
A 
A 

\ A / 



with A G C being an arbitrary constant, for the vector / e C°°(M^; C^) one obtains easily, owing 
to the relationships ()1.4[) . the following functional representation: 



fiix,t) 

f2ix,t) 



Mx,t) 



(4.10) 



gi{u -tv + wf/2 - xi^/3!, v ~ wt + zt'^/2, w - zt, 

x-tu + vt^/2 - wt^/S\ + zt'^/41), 

t\gi{u -tv + wt^/2 - xi^/3!, v - wt + zt'^/2, w - zt, 

x-tu + vt^/2 - wt^/i\ + zt'^/Al) + 

+g2{u -tv + wt'^/2 - xi^/3!, v - wt + zf/2, w - zt, 

x-tu + vt^/2 - wiV3! + zf^/Al), 

\^—gi{u -tv + wt^/2 ~ xt^jZX, v-wt + zt'^/2, w 

x-tu + vt^/2 - wt^/'i\ + zt^/4l) + 

+tXg2{u -tv + wt'^/2 - xt^/3l, v - wt + zt^/2, w - zt, 

x-tu + vt^/2 - wt^/3\ + zt^/Al) + 



zt. 



wt + zt /2, W — zt. 



fi{x,t) 



+g3{u -tv + wt72 - xt^/3\, V 
x-tu + vt^/2 - wt^ li\ + zt''/4!), 

^ ir~g\{u - tv + wt^ 12 - xt^ l'i\,v - wt + zi^/2,w- zt, 

x-tu + vt^/2 - wt^/3l + zf^/Al) + 

t^ 
+X'^—g2{u -tv + wt^/2 - xt^/?,\, v-wt + zt'^/2, w - zt, 

x--tu + vt^/2 ~ wt^/i\ + zt'^/Al) + 

+tXg3{u -tv + wt^/2 ~ xt^/il, v-wt + zt'^/2, w - zt, 

x-tu + vt^/2 ~ wt^/i\ + zt^/'i\) + 

+g4{u ~tv + wt'^/2 - xt^/3\, v-wt + zt^ /2, w - zt, 

+2 /o „,,.3 , 



x-tu + vty2 - wtysi + zr/4!), 

where gj G C°°{W^;C),j = 1,4, are arbitrary smooth complex valued functions. 



Based now on the expressions (|4.9p and (|4.10p . one can construct the related Lax type repre- 
sentation for dynamical system (|4.ip in the following compatible form: 



(4.11) 
where 



fx^e[u,v,z,w;X]f, ft:^p{e)f, p{e) := -ue[u,v,w,z;X]+q{X), 



i[u, V, w, z; A] := 



/ — A^Ua; X^Vx 

-4A'' 3X^u,j: 

-lOA^n eA-* 

V -20X^r2 lOA^ri 



~XWx 

-2}?v^ 

-iX^u^ 

-4A4 



XWx 
A^Wa; 

A3 



(4.12) 



p{e) = 



/ XuUx —}?UVx 

A -I- 4A'^M -"iX^uUx 
lOA^Mn A - QX'^u 



I 



XuWx 

2X^uVx 

iX^UUx 



q{X) := 



-UZx 
-XuWx 

-X^uv,, 



/ \ 

A 

A 

V A y 
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the mappings r^ : 7M — > M, j = 1, 2, satisfy the functional-differential equations 

(4.13) Dtri+riD^u^l, Dtr2 + r2D^u = n, 

similar to (|3.13l) . considered already above, thereby being a Lax type integrable dynamical system 
on the functional manifold Ai. 

The equations (I4.13[) . as it is easy to demonstrate [Mj [11] by means of standard differential- 
algebraic methods, possess a lot of different solutions, amongst which there are functional expres- 
sions: 

,, ,, ^ .uw'^ vw'^ vw^ 7w^ w^ , 

(4-14) n = ^.(tt-t-^t-t 



2z2 3z3 24z4 120z4 144^5' 
uw^ vw'^ 3w® """'^ ""^ 



r2 = D-A^T-T - T-r + ZWZK + 



3z3 6z4 80z5 120z5 420z6' 

Owing to the existence of the Lax type representation (|4.11l) , (|4.12p and the related gradient like 
relationship (I2.32[) . we can easily derive that the Hamiltonian system (I4.ip is also simultaneously 
a bi- Hamiltonian flow on the functional manifold A4. And as it was mentioned before, it possesses 
many bi-Hamiltonian structures, depending on a chosen solution to the corresponding functional- 
differential equations (|4.13p . 

In addition, we can construct, making use of the results above and the approach of Section 1, 
the infinite hierarchies of related conservation laws for (|4.ip . both dispersionless polynomial and 
dispersive non-polynomial ones: 



a) polynomial conservation laws: 

^27r r2-ir 



(4.15) ij(^) = / dx{vw^-uz^), 7f (1^) = / dxz^{w^ -2vz) 




H^'^^'> = / dx{3Ua;z'^ +4WxVZ + 2Za:Vw), i/^^o)^ / dx{zxW - zw^), 

Jo Jo 

p2lT 

H^ / dx[3uxz(3uz + 2vw) — 6vxz{uw + V ) + 

+6wx{uvz + 2uw'^ — v^w) + z{w^ — 2vz)], 

p2lT 

H^^^) = / dx[ki{za,{2uw-v'^) + z^) + k2((2w.x{uz-vw) + 2zx{v'^ -uw))], 
Jo 



b) non-polynomial conservation laws: 

(4.16) 

H^^°'> = / dx{wl^2vA)^ 
Jo 



(•2-K 



XX'^X ^X'^XX 1^ ^X^^XX ^XX^^XX j 5 

^ ^ 

^27r 1 

H^ ' = dx(9u^Zx ~ GuxVxWx + 2v^ — \2vxZx + &w^ j , 

^27r 
H^^^'^ = / dx\u{2VxZx - wl) + v{VxWx ~ iUxZx) + 

si 

+w{uxWx-vl + 2zx) + z{uxVx-2wx)\ ], 

rrfl4) ,_ rr(14) , ,_ rr(14) , tAU) rrflS) 7„ rr(15) , ,_ rr(15) , rr(15) 
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where 

Hf = / \Uxx{2VxZz - wl) + Vxx{v,^Wx - '?>Uj:Zj) + 

-^Wxx{ua:Wx -vl+ 2Zx) + Zxx{uxVx " 2^^;) j , 



1I2 — / uX\ZxWxx Zxx^x ] 5 



r(i4) _ 



H)^ = I dx[ki {v{2VxZx - wl) + z{AZx - UxWx) + wiVxWx - iUxZx)) + 

Jo 

+ k2z(2Zx +vl -UxWx)]'^, 



Hf = / dx[u,j.^,j.{2v^Z,j. - wl) + VxxxiVxWn, - 5u^z^) + 

Jo 

+ Zxxx{UxVx - '^Wx) + Wxxx{UxWx - vl + 2Zx) + 

(^4.1 t ) ~r^Uxx['^xxZx '^^xZxx ~r WxxWx) ~r oVxx\^'^xZxx 

-VxxWx + VxWxx) - Swl^Ux] 5 , 
^27r 1 

H2 = dxlAulwl ~ Auxvlwx — SuxZxWx + v'l + +4vlzx + 4:zl] , 

<.27r 
H^ = / dx {k3[u{z^Wxx - Z^^Wx) + 

Jo 

'y^y'^xZxx '^xxZxJ I ZZxx I WyUxxZx '^xZxx )\ 1 

-\'k4^[z['UxxWx UxWxx H~ '^Zxx) + W[UxxZx UxZxx VxxWx + VxWxx)\ + 

1 

+ hzx{vl - 2UxWx + 2z^)} ' , 

and kj G M, j — 1,5, are arbitrary constants. We observe also that the Hamihonian functional 
coincides exactly up the sign with the polynomial conservation law ij'^' e D{A4). 



,(0) 



.w 



Remark 4.1. It is worth here to remark 14 that the generalized Riemann type hydrodynamical 
equation (jl.ip can be once more naturally generalized to the following also integrable Riemann 
type equation 

(4.18) D^u = 0, Dt := d/dx + a{u)d/dt, 

where A^ e Z+ and a e C°° {M ; M) is an arbitrary smooth mapping. The corresponding to (14.181) 
nonlinear dynamical system 

= u^^^ - a{u)u^^\ 
(4.19) 

ui^-') = u(^-i)-a(^)4^-2) 
uf-^) = -a(^)4^-^^ 
will be also a bi-Hamiltonian Lax type integrable dynamical system on the phase space A4. 
Thereby, the calculations above ensue the formulation of the following proposition. 

Proposition 4.2. The Riemann type hydrodynamical system il.l\) at N = A is equivalent to 
a completely integrable bi-Hamiltonian flow on the functional manifold M, allowing the Lax type 
representation (RT771) and whose co-implectic structure is given by expression (RT^. 

Concerning the general case N G Z_|_, applying successively either the symplectic approach 
J • 1 -1 j_i_ - j'ix _L • ^ 1 ^1 1 '^ j_i 1 J • 1 •__ n~n rm _i -i__l-*__ r i j_i_ 
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the Riemann type hydrodynamical system ()l.ip and (|4.18|) the infinite hierarchies of dispersive 
and dispersionless conservation laws, co-symplectic structures and related Lax type representations, 
what is a topic of the next work under preparation. 

5. Conclusion 

As follows from the results obtained in this work, the generalized Riemann type hydrodynam- 
ical equation (jl.ll) possesses many infinite hierarchies of conservation laws, both dispersive non- 
polynomial and dispersionless polynomial. This fact can be easily explained by the fact that the 
corresponding dynamical system (II. 2[) allows many, plausibly, infinite set of algebraically indepen- 
dent compatible implectic structures, which generate via the standard gradient like relationship 
(|2.30p the related infinite hierarchies of conservation laws, and as a by-product, infinite hierar- 
chies of the associated Lax type representations. Such a situation within the theory of Lax type 
integrable nonlinear dynamical systems meets, virtually, for the first time and may appear to be 
interesting from different point of view, as well as theoretical and practical. Keeping in mind these 
and some other important aspects of the generalized Riemann type hydrodynamical equation (jl.ip . 
we consider that they deserve additional thorough investigation in the future. 

6. Acknowledgements 

Authors are sincerely appreciated to Profs. F. Calogero, M. Pavlov, M. Blaszak, Z. Peradzyhski, 
J. Slawianowski, N. Bogolubov (jr.) and D. Blackmore for useful discussions of the results obtained. 
The warm thanks go to our colleagues Dr. J. Golenia and Dr. P. Holod for instrumental help in 
editing the manuscript. The last but not least thanks go to Referees who generously mentioned 
some important points related with the integrability problem treated in the work. 



[6: 

[7: 

[9 

[lo; 
[11 
[12: 

[is: 



[14: 
[is: 

[16 

[17] 

[is: 

[19 

[20: 



References 

Whitham G.B. Linear and Nonlinear Waves . Willey-Interscience, New York, 1974, 221p 

Gurevich A.V. and Zybin K.P. Nondissipative gravitational turbulence. Sov. Phys.-JETP, 67 (1988), p.p. 1-12 
Gurevich A.V. and Zybin K.P. Large-scale structure of the Universe. Analytic theory Sov. Phys. Usp. 38 (1995), 
p.p. 687-722 

Prykarpatsky A.K., Blackmore D. and Bogolubov N.N. (jr.) Hamiltonian structure of Benney type hydrody- 
namic systems and Boltzmann-Vlasov kinetic equations on an axis and some applications to manufacturing 
science. Open systems and Information dynamics. 6, 1999, p. 335-373 

Hunter J, Saxton R. Dynamics of director fields. SIAM J. Appl. Math 51 (1991) 1498 - 521. 
Lenells J., The Hunter-Saxton equation: a geometric approach, SIAM J. Math. Anal. 40 (2008), pp. 266-277. 
BruneUi L. J. and Das A. Journal of Mathem. Phys. 2004, 45, p. 2633 
Pavlov M. The Gurevich- Zybin system. J. Phys. A; Math. Gen. 38 (2005), p. 3823-3840 
Sakovich S. On a Whitham-Type Equation. SIGMA 5 (2009), 101 

Golenia J., Pavlov M., Popowicz Z. and Prykarpatsky A. On a nonlocal Ostrovsky- Whitham type dynamical 
system, its Riemann type inhomogenious regularizations and their integrability. SIGMA 6 (2010), p. 1-13 
The differential-algebraic integrability analysis of the generalized Riemann type and Korteweg-de Vries hydro- 
dynamical equations, (in press.) 

Marsden J. and Chorin R. Mathematical backgrounds of fluid mechanics. Springer, New York, 1993 
Bogolubov N.(jr.), Prykarpatsky A., Gucwa I. and Golenia J. Analytical properties of an Ostrovsky- Whitham 
type dynamical system for a relaxing medium with spatial memory and its integrable regularization. Preprint 
ICTP-IC/2007/109, Trieste, Italy, (available at; http://pubUcations.ictp.it) 

Golenia J., Bogolubov N. (jr.), Popowicz Z., Pavlov M. and Prykarpatsky A. A new Riemann type hydrody- 
namical hierarchy and its integrability analysis. Preprint ICTP, IC/2009/095, 2009 

Prykarpatsky A.K. and Prytula M.M. The gradient- holonomic integrability analysis of a Whitham-type non- 
linear dynamical model for a relaxing medium with spatial memory. Nonlinearity 19 (2006) 2115—2122 
Prykarpatsky A.K. and Prytula M.M. The gradient-holonomic integrability analysis of a Whitham type non- 
linear dynamical model for a relaxing medium with spacial memory. Proceeding of the National Academy of 
Sciences of Ukraine, Math. Series, N.5, (2006), pp. 13-18 (in Ukrainian) 

Faddeev L.D. and Takhtadjian L.A. Hamiltonian approach in solution theory. N.Y., Springer, 1986, 476p. 
Prykarpatsky A. and Mykytyuk I. Algebraic Integrability of nonlinear dynamical systems on manifolds: classical 
and quantum aspects. Kluwer Academic Publishers, the Netherlands, 1998, 553p. 

Mitropolsky Yu., Bogolubov N. (jr.), Prykarpatsky A. and Samoylenko V. Integrable dynamical system: spectral 
and differential-geometric aspects. Kiev, "Naukova Dunka", 1987. (in Russian) 

Hentosh O., Prytula M. and Prykarpatsky A. Differential-geometric and Lie-algebraic foundations of investi- 
gating nonlinear dynamical systems on functional manifolds. The Second edition. Lviv University Publ., 2006 
(in Ukrainian) 



THE NON-POLYNOMIAL CONSERVATION LAWS AND INTEGRABILITY ANALYSIS OF GENERALIZED RIEMANN TYPE HYDRODYNAMICA 

The Institute of Theoretical Physics, University of Wroclaw, Poland 
E-mail address: ziemekOift.uni. wroc.pl 

The AGH University of Science and Technology, Krakow 30059, Poland, and the Ivan Franko State 
Pedagogical University, Drohobych, Lviv region, Ukraine 
E-mail address: pryk . anatSua . f m 



